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We propose a new method to calculate polarization induced interfacial charges in semiconductor heterostruc- 
tures using classical electrostatics applied to real-space band diagrams from first principles calculations and 
apply it to GaN/AIN heterostructures with ultrathin A1N layers (4-6 monolayers). We show that the cal- 
culated electric fields and interfacial charges are independent of the exchange-correlation functionals used 
(LDA and HSE06). We also find the calculated interfacial charge of (6.8 ±0.4) x 10 13 cm -2 to be in excellent 
agreement with experiments and the value of 6.58 x 10 13 cm -2 calculated from bulk polarization constants, 
validating the use of bulk constants even for very thin films. 



Gallium nitride (GaN) and aluminum nitride (A1N) are 
binary III/V semiconductors with wide direct band gaps 
of 3.4 and 6.2 eV, respectively, making them ideal candi- 
dates for applications in optoelectronic, high-power and 
high-frequency devices. In addition, group Ill-nitride 
heterostructures have strong spontaneous polarization 
(Psp) due to their lower-symmetry wurtzite structure 
and piezoelectric polarization (Ppe) due to their large 
lattice mismatch. This large macroscopic polarization 
leads to a high density of charge at the heterointerface. 
Polarization engineering has been successful to induce 
high sheet density electron channels without the disad- 
vantages of impurity doping, enabling new routes for de- 
vice design. Scaling of high electron mobility transistors 
for higher frequency applications has been shown to re- 
quire the use of ultrathin AlGaN and A1N cap layers with 
thicknesses between 0.5 nm and 20 nmi 2 - 

To date, polarization in heterostructures has been 
modeled using bulk constants for the calculation of the 
interface charge density and validated only for one case of 
thin GaN and A1N layers with equal thickness. 3 However, 
the validity of this approximation for arbitrary "non- 
bulk" layer thicknesses is still an open question^ 

In this letter, we propose a novel method to calculate 
interface charges, where classical electrostatics are ap- 
plied to real-space band diagrams obtained from density- 
functional theory (DFT) calculations and use it to cal- 
culate the interface charges in GaN/AIN hetero junctions 
with ultrathin A1N layers (4-6 monolayers (MLs) thick). 
We show that this method is more convenient and consis- 
tent for calculating interface charges than the commonly 
used Bader analysis. 5 Our results are finally used to ex- 
amine the validity of using bulk constants in the limit of 
small film thicknesses. 

Within the traditional bulk approximation^ the po- 
larization induced interface charges in GaN/AIN thin 
films are calculated from DFT-derived polarization 
coefficients^ and elastic constants. There, the charge den- 
sity asp induced by spontaneous polarization Psp is 



Using the values of P^ N = -0.081 cm" 2 and P|ff N = 
-0.029 cm" 2 from Ref. 7. we get a S p = 3.25 x 10 13 cm" 2 . 
Similarly, the interface charge density <tpe induced by the 
piezoelectric polarization Ppe is 

crpE = e 33 e 3 + 2e 3 iei, (2) 

with piezoelectric constants e 33 and e 3 i, in-plane and 
out-of-plane strains e\ = (ac a N — &ain)/&ain and e 3 = 
— 2eiCi 3 /c 33 , and A1N elastic constants ci 3 and c 33 . With 
Eqs. (pQ) and ([2]), the total charge density asp + <tpe is 
then independent of the layer thicknesses. 

For our DFT calculations, we used the VASP code 9 
with projector augmented wave (PAW) potentials^ 
within the local density approximation (LDA)ii as had 
been used in previous works ^Additionally, as LDA 
underestimates band-gaps of solids by 30% — 100%^ 
we repeat some band structure calculations with a 
more accurate (but computationally expensive) modi- 
fied HSE06 (mHSE06) hybrid exchange-correlation (XC) 
functional 13 for PBE 14 -relaxed structures. We use an 
optimized Hartree-Fock fraction of a = 0.35 to get good 
agreement with the experimental bandgaps of both GaN 
and A1N. A plane-wave cutoff energy of 400 (500) eV 
was used for LDA (mHSE06) calculations, along with 
T-centered Monkhorst-Pack meshes 15 with divisions Ni 
such that the product of Ni with the corresponding lat- 
tice constant ai was 25, 60, and 12 A for LDA relaxations, 
LDA band structures, and mHSE06 calculations, respec- 
tively. The Ga d-orbital electrons were treated as core 
(valence) electrons for LDA (mHSE06) calculations. 

Periodic (110) x (110) x (OOn), n = 4, ...14, superlat- 
tices along the [001] direction of the wurtzite-structure 
were used to simulate the different heterostructures. Ul- 
trathin A1N layers consisted of 4 to 6 MLs in combination 
with 4 to 24 MLs of GaN. To simulate thin films on GaN 
substrate, the in-plane lattice constants of the supercells 
were kept fixed to match that of relaxed GaN while al- 
lowing the out-of-plane constants to relax. 

The band diagram along the z-axis in real space was 
obtained by averaging over the atom-projected DOS of 
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FIG. 1. (Color online) Relaxed (22 ML GaN)/(4 ML AIN) 
heterostructure (a), its layer projected density of states (b), 
and its variation of bond lengths along the [001] direction (c). 

all atoms in the layers between z and z±2.4 A. The exact 
band edges of the projected DOS were validated with the 
help of band structure calculations in reciprocal space. 

The LDA optimized lattice parameters for GaN (a = 
3.17 A, c = 5.16 A) and AIN (a = 3.09 A, c = 4.95 
A) agree within less than 1% with the experimental val- 
ues, a = 3.16 A, c = 5.13 A, for GaN 16 and a = 3.11 
A, c = 4.98 A for AIN. 17 The bulk band gaps for GaN 
and AIN were found to be 2.13 eV and 4.53 eV, respec- 
tively, compared to experimental values of 3.4 eV and 6.2 
eV. The mHSE06 band gaps for PBE relaxed lattice con- 
stants of a = 3.22 A, c = 5.24 A for GaN and a = 3.13 A, 
c = 5.01 A for AIN are 3.61 eV and 6.04 eV, respectively. 
Our LDA values for c±s and C33 of AIN are 113 GPa and 
371 GPa, respectively, in good agreement with previous 
calculations^ and experimental values of 99 ±4 GPa and 
389 ± 10 GPaJ^ Combining them with the Ppe constants 
from Ref. 7, we find cr PE = 3.33 x 10 13 cm" 2 . Adding 
asp and ope, we find the total interface electron density 
to be 6.58 x 10 13 cm -2 , independently of the layer thick- 
nesses. In the band diagram of a heterostructure with 4 
MLs of AIN and 22 MLs of GaN (Fig. [Q, we find that 
band bending is affected by both P$p and Ppe and is 
a strong function of bond lengths. The electric field £ 
obtained from the average slope of the bands dE/dz di- 
vided by the electron charge e is found to be — 1.87 x 10 6 
V/cm and 1.11 x 10 7 V/cm in the GaN and AIN layers, 
respectively. We calculate the interface electron density 
dint from applying the integral form of Gauss' law2£ to a 
volume including the interface, resulting in 



0"int = ^GaN^G 



aN 



£ain£ 



AIN- 



(3) 



Using experimental permittivities <SQaN = 7.88 x 
10- 11 F/m and £ A in = 8.01 x lO" 11 F/m, we find 



a int (GaN/AlN) = -cr int ( AIN/ GaN) = (6.6 ± 0.3) x 10 13 
cm -2 . The error bars are due to the uncertainty in find- 
ing the slope of the bands accurately as a result of the 
very low DOS present at the band edges. From the dif- 
ferent supercells used, we find cr int to be independent of 
the thickness of the GaN layer with an average value of 
(6.7 ± 0.3) x 10 13 cm" 2 for supercells with 4 AIN MLs 
and (6.9 ±0.3) x 10 13 cm -2 for the ones with 6 AIN MLs. 
Since the two values agree within their errorbars, we take 
an average value of (6.8 ±0.4) x 10 13 cm -2 as the charge 
density at the GaN/ AIN interface, which is at the upper 
end of the experimental range of (1 to 6)xl0 13 cm -2 ^ 2 - 

In order to examine the dependence of our results 
on the value of the band-gap, we repeat the interface 
charge density calculations for (4 ML AIN/ 12 ML GaN) 
and (6 ML AIN)/ (8 ML GaN) heterostructures with the 
mHSE06 functional. We find mHSE06 charge densities 
of (6.7±0.3) x 10 13 cm" 2 and (6.6±0.3) x 10 13 cm" 2 , re- 
spectively, as compared to LDA values of (6.4±0.3) x 10 13 
cm -2 and (6.8 ±0.3) x 10 13 cm -2 . Thus, the influence 
of the XC functional used on the slope of the bands in 
these heterostructures is found to be negligible within 
the computational error bars, allowing the calculation of 
interface charges even for relatively large systems (with 
100s of atoms) using the much faster LDA functional. 

In order to corroborate our calculations, we determine 
the variation of the electric field with GaN layer with an 
alternative approach. The periodic boundary conditions 
in our calculations along with the conservation of charges 
require that the integrated electric field, or overall elec- 
trostatic potential, is zero. Thus, Eq. (j3j) requires that 
the charges at the GaN/ AIN and AIN/GaN interfaces are 
of equal size, but opposite sign. Since the average slopes 
of the band edges in the GaN and AIN layers away from 
the interface are constant, this condition then leads to 



^GaN^GaN = ^AIN^AIN, 



(4) 



where fcaN and £ain are the thicknesses of the GaN and 
AIN layers, respectively, measured from one of the inter- 
faces. Combining Eqs. (j3j) and (j4]), we find 



^GaN = ^int/ [£GaN(l ± ^GaN^A^)] • 



(5) 



Figure [2] shows a plot of the electric field in the GaN 
layer vs. tc a N and a corresponding fit of Eq. (J5j) for het- 
erostructures with 4 and 6 AIN MLs. The error bars 
are due to the uncertainty in delineating the band edges. 
With £ain fixed, the only adjustable parameter in Eq. (0 
is dint. Equation (j5]) fits the calculated electric field val- 
ues well with a least-squares value of <Ji nt = 5.85 x 10 13 
cm -2 . By extrapolating the curve to tc a N = 0, we get 
the limiting value of the electric field in the GaN layer to 
be 1.19 x 10 7 V/cm (with a negligible difference of 0.03 
MV/cm between the supercells with 6 and 4 MLs of AIN). 
Both values are in reasonable agreement with the values 
determined above from Eq. J3},(6.8 ± 0.4) x 10 13 cm -2 
and (1.38 ±0.08) x 10 T V/cm, respectively, thus demon- 
strating the consistency of our approach. 
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FIG. 2. (Color online) Variation of electric field in the GaN 
layer with changing GaN layer thickness for two fixed A1N 
interlayer thicknesses of 4 (blue squares) and 6 (red circles) 
MLs from LDA calculations. The lines are fits using Eq. (5). 



As a further check, we also calculate the interface 
charge density directly from the DFT charge distribution 
at the interface, 3 For that, we use the Bader approach, 21 
where the number of electrons on each atom is deter- 
mined by integrating the charge density over the volume 
enclosed by the atom's zero flux surface. Charges on 
atoms at defects can then be determined relative to the 
charge on perfect bulk atoms, which is 1.85e and 1.06e 
for N atoms in bulk A1N and GaN, respectively. Refer- 
ence charges for interfacial N with 3 Al and 1 Ga bonds 
and vice versa are calculated from a structure with alter- 
nating GaN and A1N MLs with values of 1.61 e and 1.1 9e 
respectively (due to the lack of any other conventional 
definition). This treatment gives the average interface 
charge density to be 5.95 x 10 13 cm -2 and 6.50 x 10 13 
cm -2 for the heterostructures with 4 MLs and 6 MLs 
of A1N (with different monolayers of GaN), respectively 
which is close to the values that we got from the electro- 
statics treatment. Again, since there is not a significant 
change in the charge density with 4 and 6 MLs of A1N, we 
take an average value of 6.23 x 10 13 cm -2 as the interface 
charge density in GaN/AIN heterostructures. 

One of the central questions of this letter concerns the 
validity of using bulk polarization constants to calculate 
interface charges in the limit of ultrathin layers. Since 
the value of interface charge density for the GaN/AIN 
heterostructures with ultrathin A1N layers obtained with 
our approach, (6.8 ± 0.4) x 10 13 cm -2 , agrees well with 
the value of 6.58 x 10 13 cm -2 from bulk polarization 
constants, our results suggest that the use of bulk polar- 
ization constants to calculate interface charges is a valid 
approach even for heterostructures with very few MLs. 

In conclusion, we have calculated the interfacial po- 



larization charge density in GaN-AIN heterostructures 
from a combination of DFT with electrostatics and find 
0"int = (6.8 ± 0.4) x 10 13 cm -2 , at the upper end of the 
experimental range of (1 to 6)xl0 13 cm~ 2 i^ This value 
is in the same range as the charge densities determined 
from the dependence of the electric field in GaN on the 
layer thickness, 5.85 x 10 13 cm -2 , and from a Bader anal- 
ysis, 6.23 x 10 13 cm -2 . The latter however has uncer- 
tainties about interfacial reference charges and number 
of layers to include into the charge density calculation. 
In comparison, we find that the combination of ab-initio 
real-space band structures and electrostatics eliminates 
these uncertainties and is applicable to a wide range of 
systems beyond the AIN/GaN system studied. Our value 
agrees well with the interface charge density calculated 
using bulk constants of 6.58 x 10 13 cm -2 , thus validat- 
ing the use of the former in heterostructures where layer 
thicknesses approach "non-bulk" dimensions. Addition- 
ally, this method could also be used to calculate charges 
at more complex interfaces, such as, among others, in 
digital superlattices and tunnel barrier systems. 
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